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ABSTRACT 

We present photometry and moderate-resolution spectroscopy of the luminous red 
variable [HBS2006] 40671 originally detected as a possible nova in the galaxy M33. 
We found that the star is a pulsating Mira-type variable with a long period of 665 
days and an amplitude exceeding 7 mag in R band. [HBS2006] 40671 is the first 
confirmed Mira-type star in M33. It is one of the most luminous Mira-type variables. 
In the K band its mean absolute magnitude is Mi<- = —9.5, its bolometric magnitude 
measured in the maximum light is also extreme, Mf,,,; — —7.4. The spectral type of 
the star in the maximum is M2e — M3e. The heliocentric radial velocity of the star is 
—475 km/s. There is a big negative excess (—210 km/s) in radial velocity of [IIBS2006] 
40671 relative to the average radial velocity of stars in its neighborhood pointing at 
an exceptional peculiar motion of the star. All the extreme properties of the new Mira 
star make it important for further studies. 

Key viTords: galaxies - optical: variable stars. 



1 INTRODUCTION 

The General Catal ogue of Variable Stars (GOVS) 
l|Kholopov at aL 19851 ) gives a determination of Mira (Omi- 
cron) Ceti-type variables: "These are long-period variable 
giants with characteristic late-type emission spectra (Me, 
Ce, Se) and variation amplitudes from 2.5 to 11 magnitudes 
in V . Their periodicity is well pronounced, and the periods 
lie in the range between 80 and 1000 days. Infrared ampli- 
tudes are usually less than those in the visible and may be 
< 2.5 mag. For example, in the K band they usually do not 
exceed 0.9 mag." Due to the large amplitude of variability 
exceeding 2.5 mag in the V band, Mira variables in nearby 
galaxies may be confused with optical novae. The confusion 
can be solved, if subsequent spectral or photometric inves- 
tigations reveal a Mira-type star with a cool M-type spec- 
trum and a long periodic variability. Recently, peculiar red 
novae were introduced as a new class o f astrophysical ob- 
jects (e.g., iGoranskii fc Barsukova 200?! ). It is represented 
by stars such as V838 Mon, V4332 Sgr and V1006/7 in M 31. 
Nova Sgr 1943 (VI 148 Sgr) may possibly also belong to this 
class. Their remnants may be cool luminous L- or M-type 
supergiants, which in the course of discovery may be con- 
fused with both Mira-type stars and classical novae. Nova 
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Sgr 1943 reached a majci mum photom etric brightness of 8.0 
mag, it was described bv lMavall (19491 ) as a late K-type star 
with TiO bands in the spectrum. It was found that the star 
is not of Mira-type, and later the star was lost. 

In a search for novae in nearby galaxies one may run 
into similar problems. A cool supergiant may be misidenti- 
fied as a nova when an unfiltered CCD image is compared 
with B ox V images. The M31 Nova candidate 2008-09b 
l|Barsukova et al. 20081 ) turned out to be a star with M9- 
type spectrum. It may be a red nova or even a red su- 
pergiant showing no variability at all. There are two addi- 
tional nova candidates in M31 that t urned out to be lum i- 
nous Mira-type st ars: M31N 19 95-lle dShafter et al. 2008al ). 
M31N 2007-llg (|Shafter et al. 2008b ). The Mira-type star 
in IC 1613 described by Kurtev et al. ( 2001*) has been dis- 
covered as a nova as well (jKing et al. 1999 ). In this paper 
we study a Mira-type star in the galaxy M33, which has 
been discovered as a nova candidate. 

An outburst of a source with coordinates RA = 
0l''34"'27^ 13, Dec = 30°58'42'.'7 (J2000) in M33 was de- 
tected by iNishivama fc Kabashima (20091 ). They reported 
it as a possible nova (with 18.6 mag) detected in five frames 
taken around 2009 August 17.807 UT exposed by 60 sec 
using the 0.40-m f/9.8 reflector and unfiltered CCD of the 
Miyaki Argenteus Observatory, Japan, reaching a limiting 
magnitude of 20.2. Nothing was visible at this location on 
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Figure 2. Periodogram of [HBS2006] 40671 calculated in 
the period range bet ween 200 and 7000 days using the 
iLafler fc Kinman fl965h phase dispersion minimization method. 
6 is a normalized dispersion parameter. 



Figure 1. Finding chart of [HBS2006] 40671 made from the R 
band Subaru telescope frame taken on 2002 November 4. The 
new Mira-type star is marked by two lines. Size of the image is 
25" X 25". 



their previous frames taken on 2008 December 25.559 UT 
and 2009 August 07.737 UT with limiting magnitudes of 
20.0 and 19.4, respectively. The source was confirmed by 
these authors on 2009 August 18.696. But later they with- 
draw the interpretation as possible nova as it did not change 
brightness. 

The source was identified by us with a star number 
40671 in the photometric survey of variable stars in M 33 
carried out bv llfartman et al. f2006l ). We will name it here- 
after [HSB2006] 40671. The survey observations were carried 
out between 2003 August and 2005 January. The authors 
classify the star as a Long Period Variable (LPV). During 
their observations (filled circles in Fig. 3), the star was red 
with R — I ^ 2.7 and ~ 4.1 mag at maximum and min- 
imum brightness, respectively. The star became extremely 
red during the registered minimum (R > 25 mag) . 
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Figure 3. Light curve of [HBS2006] 40671 folded modulo the 
665-day period in R ba nd. Filled circles are the observations 
by Hart man et al. (20061) . Big open squares indicate observa- 
tions from Calar Alto, the asterisk from the Subaru telescope, 
open triangle from SAO BTA, filled triangle from SAO 1-m tele- 
scope, small open squares are Sloan r' magnitudes corrected 
by 0.6 mag. Open rhombs are unfiltered observations taken by 
llNishivama fc Kabashima 20091) (error bars ~ 0.2 mag are not 
shown for clarity). At R > 25 mag the accuracy of the estimates 
is not better than 0.4 mag. 



2 OBSERVATIONS 

We have carried out spectral observations with the Rus- 
sian 6-m BTA telescope usi ng the SCORPIO spectrograph 
ijAfanasiev fc Moiseev 20051 ). The spectrum has been taken 
on 2009 October 9.874 UT in a range 4300 - 7880 A with 
spectral resolution of 13 A. The spectroscopy was followed 
by BVRc images on the same date. We found the star at 
a level of 18.8 mag in the Rc band. One more photometric 
BVRc observation has been taken with the 1-m Zeiss tele- 
scope of the Special Astrophysical Observatory (RAS) on 
2009 October 20.9 UT. 

We have searched for additional images of this star in 
Internet archives of different observatories. Dozens of images 
were found, most of them come from medium-size telescopes 
equipped with wide field mosaic CCD cameras. The star po- 
sition is well placed on the frames of the Calar Alto Obser- 



vatory taken in 1989 and 1998 with the 80/120-cm Schmidt 
telescope, on the frames taken in 2001 and 2002 in BVRI 
bands with the 8.2-m SUBARU telescope of the National As- 
tronomical Observatory of Japan. One V band observation 
from the 4.2-m William Herschel Telescope, and three SDSS 
r' band observations from the 2.54-m Isaac Newton Tele- 
scope were used as well, both telescopes belong to the Isaac 
Newton Group of Telescopes at the Roque de Los Muchachos 
Observatory at La Palma, Spain. We also apply one / band 
measurement obtained from the 4.0 -m Mayall telescope im- 
age taken in 2001 and published bv lMassev et al. f 20061 ). 

Archival images from the 8.2-m SUBARU telescope, 
the 4.2-m William Herschel Telescope, and the 2.54-m Isaac 
Newton Telescope were downloaded as raw FITS files and 
then processed in the same way as our own 6-m BTA and 
1-m Zeiss telescope images. Standard reduction procedures 
for raw CCD images were applied (bias and dark-frame sub- 
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Table 1. Photometric Observations of [HBS2006] 40671. 
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For the plot of the hght curve (Fig. 3), we subtracted 0.6 mag from Sloan r ' 
magnitudes to tra nsform them into R magnitudes (using color transformations from 
iJordi et al. (20061) and assuming an average color index V — R = 1.2 mag for the 
Mira-type star). 



tract and flat-field correction) using the SIM^ program. 
Reduced images of the same series were co-added to im- 
prove the S /N ratio and then used for photometry. We used 
GAIA0 to perform "Optimal photometry" (based on fitting 
of PSF profiles). B,V,R, and / magnitudes for compari- 
son stars were taken from the Lo cal Group Galaxy Survey 
(LGGS) catalog of stars in M33 (jMassev et al. 20061 ). 

In the case of SDSS r' band images, we computed r' 
magnitude s for comparison sta rs from BVRI magnitudes 
taken from lMassev et al. (2006l l using empirical color trans- 
formations between the SDSS u' q'r'i' z' system and Johns on- 
Cou sins U BVRI system published b y IJordi et al. (20061 ). 

iNishivama fc Kabashima (20091 ) continued their unfil- 
tered observations at the Miyaki Argenteus Observatory, 
Japan, from 2009 August 18 to 2010 January 3 and kindly 
provided us with 16 new images taken in this period. All the 

^ http : //ccd.mll . cz/| 

^ http : //www. starllnk.rl . ac .uk/gaia| 



photometric measur ements of [HBS20061 4 0671 (collected in 
addition to those of iHartman et al. 20061 ') are presented in 
Table 1. 



We created a finding chart for the star from high- 
quality R band image taken with the Subaru telescope 
on 2002 November 4 (Fig. 1). A relatively bright visual 
companion is locate d 2.1 arcsec north-east of the variable. 
iMassev et al. (20061 ) cataloged the position of this star as 
RA = 01*'34"™2r:l9, Dec = 30°58'44'.'4 (J2000) and give 
B, V, R, and / band magnitudes of 21.40, 21.11, 20.94, 
and 20.71, resp ectively. The observations of the LGGS 
(jMassev et al. 20 06) were obtained from October 2000 to 
September 2001. The variable star is not included in the 
catalog, as it is only detected on the I band image taken 
during the survey on 2001 September 18 (criterion for inclu- 
sion in catalog is detection in the three bands B, V, R). 
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Figure 4. The spectrum of [HBS2006] 40671 taken with BTA on 2009 October 9 near the maximum light. 



3 THE PERIOD SEARCH AND THE LIGHT 
CURVE 

Most of the data in Table 1 relate to the light curve 
maximum. Both ascending and descending branches of 
the light curve are not fully covered. The observations 
during these phases come mostly from the survey by 
iHartman et al. (20061 ). The survey presents observations 
when the star had a brightness between 22 and 25 
mag in the R band. The accuracy of these observations 
l|Hartman et al. 20061 ) varies from 0.04 mag at i? = 22 mag 
to 0.37 mag at 7? = 25 mag. 

To search for periodic ities we used a method proposed 
bv lLafler fc Kinman which is based on the phase dis- 

persion minimization. In this search we used only the obser- 
vations in R band because they represent the best sampling 
among the total data set. 

The periodogram was calculated in the period range 
between 200 and 7000 days (Fig. 2). The output parame- 
ter there is the inv erse dispersion value 1/6 described by 
iLafler fc Kinman (19 65). maxima in the periodogram corre- 
spond to dispersion minima. There are essentially no peaks 
with period values less than 200 days. This range of peri- 
ods is not shown in the figure. A period of 665 days shows 
the highest amplitude. There are two lower amplitude peaks 
at 3500 and 406 days (Fig. 2). However, these periods give 
much worse light curves. Longer periods contradict with 
the fast changes of brightness on ascending and de scend- 
ing branches of the light curve (|Hartman et al. 20061 ). The 
phased light curve in the R band is presented in Fig. 3. It 
was calculated with the following light elements 



Max = 2455080 + 665 * E. 



(1) 



Figure 3 shows that we do not resolve and possibly not even 
cover the minimum of the light curve. The star in minimum 
brightness is not brighter than R = 25.4 mag. The total 
amplitude of the variability is bigger than 7 mag in the R 
band. Also the ascending branch is not well covered. There- 
fore, the epoch of the maximum can only be constrained to 
±5 days. However, it is obvious that the light curve of this 
Mira variable is asymmetric, the ascending branch is steeper 
than the descending one. 

We tried to refine the maximum in the light curve us- 



ing archival observat ions of the infrared satellite Spitzer 
uinn et al. 20071 ). In five observations on 2004 Jan- 
uary 09, July 22, August 16, and 2005 January 21 and Au- 
gust 25, the star varied in the range of 14.24 - 14.97 mag, 
14.16 - 14.64 mag, and 13.66 - 14.04 mag in the 3.6, 4.5, 
and 8 \im bands, respectively. A maximum is found during 
July /August 2004 corresponding to phases in the interval 
0.18-0.22 according to Equation 1. Therefore the infrared 
maximum seems not to be in phase with the optical maxi- 
mum exactly, but it is very close to it (see Fig. 3). 



4 SPECTRUM 

The calibrated spectrum of [HBS2006] 40671 in the 4300 
- 7880 A spectral region is shown in Fig. 4. Molecular 
TiO bands dominate in the spectrum. We clas sify the star 
as an oxygen-rich (M-t ype or O Mira) star (|Smak 19641 : 
Ijaschek fc Jaschek 19871 ). Strong narrow emission lines of 
hydrogen are present as well; this is typical for Mira stars 
during maximum light. Using the spectrum we estimate 
a contribution of the TiO bands into the star brightness 
in B and V bands. They are 5{B) — —0.28 mag and 
5{V) = —0.16 mag. These values agree with the spectral 
type Ml - M2 jSmak 19641 ). 

We have also determined the spectral type of the 
star using [TiO ji and [TiO] 2 indices introduced by 
lO'Connell (197^ ) and measured in our calibrated spectrum. 
The spectral indices were derived from fiuxes measured in 
30 A bandpasses centered at three wavelengths each for 
[TiO]i (6125, 6180 and 6370 A) and for [TiO] 2 (7025, 7100 
and 7400 A) as [TiO]i = -HO. 52 and [TiO] 2 = +0.55, respec- 
tively. These index values correspond to spectral class M3. 
Such a spectrum is natural for Mira variables in maximum 
light, while in the minimum they have usually later spectra. 

The Ha emission hue of [HBS2006] 40671 is very nar- 
row and not resolved in our spectrum. Its equivalent width 
is EW = 9 A. For Ha we measure a heliocentric radial ve- 
locity of —440 ± 15 km/s. That this high radial velocity 
is not caused by calibration problems was verified by mea- 
suring the radial velocity of the [OI] A6300 sky line which 
was determined as km/s as expected. Also, at a distance 
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of - 1/6 from [HBS200 6] 40671, by chance the high ion- 
ization HII region Z378 (|Courtes et al. 1987^ was captured 
in the slit. The heliocentric radial velocity of its Hq line 
is —280 ± 10 km/s. This velocity is in good agreement with 
that of the HI radio emi ssion at the position (— 265±5 km/s, 
iReakes fc Newton 197^ . The difference of 1.6' between the 
HII region and the star's location in the galaxy body may 
produce a radial velocity shift of 5 - 10 km/s. Therefore 
we find a strong negative excess in the Ha radial velocity 
of [HBS2006] 40671 of -175 km/s over the galaxy radial 
velocity. 

To check whether this negative excess of the Ha ra- 
dial velocity also is reflected in other features of the spec- 
trum of [HBS2006] 40671, we cross-correlated it with spec- 
tra of known standard stars. We used two red variables, 
V934 Her (HD 154791, 4U 1700-^24, M2 HI) and KK Per 
(HD 13136, M2 lab-Ib). The cross-correlation region was 
5850 - 6700 A excluding H^ emission and interstellar Nal 
lines. The spectrum of V934 Her was taken with the SAO 
1-m Zeiss telescope o n 2006 April 6.95 UT. According to 
iGaUowav et al. (20021 ) it is an X-ray binary with a period of 
P = 404 ±3 days. The radial velocity amplitude of the pri- 
mary is only K = 0.75 ±0.12 km/s. The systemic velocity of 
V934 Her is — 48.7±0.1 km/s. The cross-correlation resulted 
in a radial velocity of —474 km/s for the Mira star. We ob- 
tained almost the same value (—477 km/s) using the second 
standard star KK P er (its radial velocity is — 39.4±0.4 km/s, 
iMarrese et al. 2003^ . 

We conclude that the star has a big negative peculiar 
radial velocity compared to the body of the galaxy (HI data) 
of —210 ± 10 km/s. Its Hq emission is redshifted relative to 
the star by -1-35 ± 10 km/s. Such a shift might be explained 
by stellar pulsations. However, it could also be caused by the 
superposition of the emission line on an absorption line that 
is not resolved in our spectrum. The huge peculiar velocity 
can not be explained by the star binarity. To get an orbital 
velocity of ~ 200 km/s in a few solar mass binary the orbital 
period must be less than 10 days, and the orbital separation 
less than 0.2 a.u. This is impossible due to the huge sizes of 
Miras. 

We tried to investigate if [HBS2006] 40671 is positioned 
in front of or behind the disk of M 33 which would indicate 
motion away from or to the disk, respectively. However, the 
signal to noise of the spectr um in the range of expected 
diffuse interstellar bands (e.g. ICordiner et al. 2010l ) was not 
sufficient. 



5 DISCUSSION AND CONCLUSION 

According to our data from the 6-m and 1-m telescopes, 
the observed photometric colours of [HBS2006] 40671 during 
maximum are B — V = 1.6 and V — R — 1.2 mag. Taking into 
account the distance modulus of M 33 (m — M)o = 24.92 ± 
0.12 mag (the distance 964 ± 54 kpc. iBonanos et al. 200^ ) 
and the Galactic extinction v alue Ay = 0.15 ma g {Ak = 
0.02 mag) in that direction (|Schlegel et al. 19981 ) we find 
an absolute magnitude of the star near maximum bright- 
ness (V= 20.02 ± 0.05 mag at the pulsation phase 0.05) 
of Mv = —5.05 ± 0.15. This is consistent with luminos- 
ity classes Ib-Iab. For an M3-typ e star the bo lometric cor- 
rection amounts to —2.34 mag (|Smak 19661 ). This leads 



to a bolometric magnitude in the brightness maximum of 
Mfco! ~ —7.4. Unfortunately, we have not enough data to 
estimate V band pulsation magnitude or mean magnitude. 

2MASS data (Skrutskic C!t_al._ 2006) show [HBS2006] 
40671 on 1997 December 5, in the pulsation phase of 0.54, 
at a brightness of 17.30 ±0.23, 15.89 ±0.16 and 15.45 ±0.15 
mag in J, H and K, respectively. There are more obser- 
vations of the star in the JHK bands taken from 2005 
September 29 to December 16 (pulsation phases 0.84-0.96) 
with the 3.8-m U KIRT telescope equipped with WFCAM 
(jCioni et al. 200^ ). The brightness of the star during this 
period was 16.86, 16.24 and 15.66 mag in the JHK bands, 
respectively. These data confirm the typical Mira behaviour, 
i.e. the decrease of pulsation amplitude with increasing 
wavelength. However we have not enough data to study 
the infrared light curve in detail. iKanbur et al. fl997l ) found 
that the K band amplitudes of O Mira variables is smaller 
for shorter periods. From the relations in that paper, the 
O Mira variables with periods longer than 420 days should 
show total K amplitudes smaller than 0.4 mag. We use the 
2MASS data for a luminosity estimate. We find an absolute 
magnitude of the [HBS2006] 40671 as Mk ~ -9.5 with a 
possible error of ± 0.15. A pplying a bolometric co rrection 
of +3.3 mag determined bv lBessell fc Wood (1983 ) to this 
K b and magnitude we find Mboi = —6.2. 

iFeast et al. (198*^ ) present period - luminosity (PL) re- 
lations for LMC Mira variables, where they found that the O 
Mira variables wit h periods longer than ^ 420 days are over- 
luminous (see also lHughes fc Wood 1990l ). Fitting the over- 
luminous branch of the long-period Mira variables (without 
the star C2 whose status is less certain, iFeast et al. 198i^ ) 
we find for our star with the 665-day period that the ex- 
pected mean luminosity is Mk ~ —9.47 mag. This agrees 
well with our estimate from the 2MASS data. Using PL 
relations for bolometric luminositie s for the overlumi nous 
long-period Mira variab les both from lFeast et al. (19891 ) and 
iHughes fc Wood (199d ) we find very similar results, the 
mean expected luminosity for the 665-day star is Mtoi ~ 
—6.5 mag. This value is close to Mboi = —6.2, found above 
from the 2MASS data taken in the minimum. On the basis 
of our spectral classification, we determined the bolometric 
luminosity of « —7.4 mag is in the brightness maximum. 
The bolometric luminosities of 0-type Mira variables are 
notably changin g during the pulsatio n cycle. Using the re- 
lations found bv lKanbur et al. (19971 ) one may expect that 
in O Mira variables with periods longer than 420 days the 
difference between maximal and mean bolometric luminos- 
ity is larger than 0.5 mag. We conclude that luminosities of 
[HBS2006] 40671 found by us are in line with the statistical 
relations, but the star's parameters put it in the extreme 
tail of these relations. 

In the presence of the huge amplitude of [HBS2006] 
40671 in the R band, not less than 7 mag, one 
may expect an even larger amplitud e in the V band 
(e.g. iBarthes. Chenevez. fc Mattel 1996t ). In the long-period 
Mira variables the amplitude dispersio n is very high 
(|Mattei et al. 19971 : [Hughes fc Wood 199d ). A visual ampli- 
tude in a Mira sta r with a period P ~ 665 days may be as 
large as 8. 3 mag (iMattei et al. 19971 ). Using relations pre- 
sented by iHughes fc Wood (1990l) we estimate a mass of 
this star as M ~ 4 solar masses. Based on its luminosity 
and the pulsation period, [HBS2006] 40671 is similar to the 
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Mira-type variable discovered bv lKurtev et al. (20011 ) in the 
IC 1613 galaxy (Mk = -9.62 mag and P = 640.7 day). How- 
ever the amplitude of the IC 1613 Mira is only 2.5 — 3 mag 
in R band. In the PL plane for Mk luminosities for LMC 
the star [HBS 2006] 40671 is locat ed in the same place as the 
IC 1613 Mira (|Kurtev et al. 20011 ). Both stars are located in 
the zone of the first-overtone pulsating Mira variables. 



IWood et al. (19831 ') found that long-period variables are 
grouped into two classes: core- helium- burning supergiants 
which are brightest, and AGE stars which are at least 1 mag 
fainter at a given period. The supergiant LPVs form a dis- 
tinct PL relation and they have lower amplitudes. The PL 
relation for si x core-helium-burnin g supergiants in M 33 was 
presented by iKinman et al. (19871 ) in their Fig. 8. The K 
magnitude for supergiants with the 665 day period is ex- 
pected to be 13.9. The 2MASS K magnitude of [HBS2006] 
40671 is 1.5 mag fainter, what suggests that this object is 
an AGB s tar. AGE stars in M 33 and M31 are discu ssed in 
papers bv ljavadi et al. (2010l ) and lRich et al. (19931 ). 



IWhitelock (20101 ) discussed a type of long-period vari- 
ables which are undergoing a hot bottom burning (HBB) in 
the base of their convective envelopes. They lie above the 
PL relation of AGB Mira-type variables. In LMC they are 
large amplitude variable stars additionally showing a lithium 
and s-process enrichment. As an e xample, she refers t o 641- 
day Mira in IC 1613 studied bv iKurtev et al. (200il ) that 
locates above the PL relation. [HBS2006] 40671 also locates 
0.35 mag above the PL rela tion for AGB stars wit h peri- 
ods over 400 days derived bv lHughes fc Wood (1990) along 
with three HBB stars studied bv lWhitelock et al. (2003 ) in 
LMC. We did not find lithium or other s-process elements 
in our low-resolution spectrum, but the luminosity excess is 
evident as that in the IC 1613 Mira variable. 

We conclude that the star [HBS2006] 40671 in M 33 is 
O Mira-type variable with extreme properties. Its pulsation 
amplitude is not less than 7 magnitudes in R band, the 
period of ^ 665 days is one of the longest known for Mira- 
type variables. The mean absolute magnitude of the star 
in K band is }Ak ~ —9.5. In maximum light its bolometric 
magnitude is estimated as M-boi ~ —7.4. It shows a spectrum 
of type M2e-M3e. 

There is a strong negative excess of —210 km/s in the 
star velocity relative to the radial velocity of the star's loca- 
tion projected to the gal actic disk. T his is a peculiar motion 
of the star. According to lFeast (20071 ) the velocity dispersion 
of Galactic 0-Miras is in the range 30 to 80 km/s. They be- 
long to the galactic populations ranging from the thin disk 
to the extended disk. In spite of the enhanced velocity dis- 
persion of the Galactic Mira-stars the peculiar velocity of 
~ -200 km/s in the [HBS2006] 40671 is rather big. 

The radial velocity of the narrow Hq emission line is 
-1-35 ± 10 km/s relative to that of the star itself. Such a 
shift in the Ha emission is typical for Mira stars. It can 
be explained by stellar pulsations and shocks in an expand- 
ing atmosphere. All the properties make the new Mira star 
important for further studies. [HBS2006] 40671 is the first 
spectroscopically confirmed Mira star in M33. 
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